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ABSTRAGT 

The observed properties of high redshift galaxies depend on the underlying foreground distribution 
of large scale structure, which distorts their intrinsic properties via gravitational leasing. We focus 
on the regime where the dominant contribution originates from a single lens and examine the statis¬ 
tics of gravitational leasing by a population of virialized and non-virialized structures using sub-mm 
galaxies at z ^ 2.6 and Lyman-break galaxies at redshifts z ^ 6 — 15 as the background sources. We 
quantify the effect of leasing on the luminosity function of the high redshift sources, focusing on the 
intermediate and small magnifications, /i < 2, which affect the majority of the background galaxies, 
and comparing to the case of strong leasing. We show that, depending on the intrinsic properties of 
the background galaxies, gravitational leasing can significantly affect the observed luminosity function 
even when no obvious strong lenses are present. Finally, we find that in the case of the Lyman-break 
galaxies it is important to account for the surface brightness profiles of both the foreground and the 
background galaxies when computing the lensing statistics, which introduces a selection criterion for 
the background galaxies that can actually be observed. Not taking this criterion into account leads to 
an overestimation of the number densities of very bright galaxies by nearly two orders of magnitude. 
Subject headings: lensing 


1. INTRODUCTION 

Detection of the high-redshift galaxies is a primary 
frontier in observational cosmology. Sampling and an¬ 
alyzing the properties of different types of high-redshift 
sources will constrain galaxy formation and star for¬ 
mation histories at different epochs, and explain their 
role in reionization and metal enrichment o f the Uni¬ 
verse (iLoeb fc Furlanettol l2013t lEllisI 120141 . For in¬ 
stance, dusty star-forming galaxies at a redshift range 
z ^ 2 — 4 are the most luminous galaxies at that 
epoch and host a considerable frac t ion of star forma¬ 
tion at z > 2 (iHughes et al. 1 119981 : iBlain et al. 1 119991 : 
iGhaoman et al. 1 1200511 . which includes t he epoch when 
the s tar-formation rate density peaked (iMadau et al. I 
Il996fl . Despite their leading role in the history of galaxy 
formation, this population of galaxies has not been well 
studied yet due to dust obscuration, and is surveyed at 
present at sub-millimeter wavelengths. A reliable num¬ 
ber counts of resolved sources from this population are 
being provided for the first time (|Karim et al. 1 1201311 
by the Atacama Large Millimeter/submillimeter Array 
(ALMA). Another population of high redshift sources 
which is not well constrained at the moment are the 
galaxies existing during the epoch of reionization and ob¬ 
served in their rest-frame UV at redshifts out to z ^ 10 
by the Wide Field Ga mera 3 infrared channel (WFG3/IR, 
iKimble et al. I (|2008[ 11 on board the Hubble Space Tele¬ 
scope (HST), with plans to push this frontier to even 
higher redshifts with the James Webb Space Telescope 
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(JWST). The identification of these galaxies is done 
based on the Lyman-break technique which relies on the 
absorption of ultra-violet photons at wavelengths short- 
ward of the redshifted Ly-a line due to neutral hydro¬ 
gen fraction. Because the UV luminosity of the Lyman- 
break gala xies strongly correlates with the sta r forma¬ 
tion rates (|Madau. Pozzetti fc Dickinson IIT998I 1. estab¬ 
lishing the UV luminosity function at high redshift is an 
essential step towards measuring the halo abundances 


at these redshifts (e.g., via a b undance-matching tech¬ 
niques (i Kravtso v et al I 12 004 iT asitsiomi et al. I 12004 
I Vale fc Ostriker 112004 iGonrov et al.ll2006fB and assess¬ 
ing the impact of these galaxies on the reionization of 
cosmic hydrogen in the first billion years after the Big 
Bang. 

The observable properties of any high-redshift popula¬ 
tion of sources, such as their abundance and luminosity 
function, differ from the intrinsic ones, since their radia¬ 
tion is subject to gravitational lensing. Along a random 
line of sight, photons are deflected on their path from 
the distant galaxies to th e observer, a process which can 
be described statistically (iTurner et al. III9^lPei II1993L 


I995I: Perrotta et al. Il200a Negrell^ 2007 : Lima et al. I 


2010a: Lima. Jain fc Devlin Il2010bl: I.Tain fc Lima 1120111: 


Ward^ 


ow et al.l 120131 1. In rare cases when the minimal 


separation between the center of the foreground struc¬ 
ture and the trajectory of the photon (referred to as its 
impact parameter) is small enough, gravitational lensing 
has a strong effect on the light rays leading to a signif¬ 
icant magnification of the background galaxy, multiple 
images and a major distortion of the source appearance 
on the sky. The photon deflection has two main impacts 
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on the observed sources: the fluxes are magnified and the 
solid angle within which the sources are observed is in¬ 
creased. Both effects need to be taken into account when 
inferring the intrinsic properties of a background popu¬ 
lation from the observed ones. The efficiency of lensing 
depends on the properties of both the foreground and 
the background populations. For instance, for the same 
set of lenses the main contribution to strong lensing with 
magnification, above 2 originates from structures at a 
mean redshift of z ~ 0.6 when the sources are the sub¬ 
millimeter galaxies (SMGs) located at z ~ 2; while the 
Lyman-break galaxies (LBGs) at z ^ 6 — 10 are most af¬ 
fected by foreground objects at z ~ 1. However, for any 
background population the strongest effect of lensing is 
expected to be on the rare brightest sources. This is 
because an observed bright source has an enhanced like¬ 
lihood of being a magnified intrinsically fainter galaxy, 
which are much more abundant, thus leading to an over¬ 
estimation of the intrinsically bright population and dis¬ 
torted appearance of its luminosity function. 

Lensing appears to be a useful tool when count¬ 
ing the numbers of high redshift galaxies and mea¬ 
suring their luminosities. The magnification makes 
the selection of strongly lensed galaxies easy for large 
area galaxy surveys. For instance, more than 85% of 
dusty sub-millimeter galaxies observed by the Herschel 
Space Observatory were confirmed to be gravitationally 
lensed by an intervening foreground structure along the 
line of sight (iNegrello et al. 1120101 : iWardlow et al.l[2(H^ 
iBussmann et al.T 201.81 1. hr addition, faint galaxies be¬ 
hind massive foreground structures can be magnified 
above the detection limit. When background samples 
behind known foreground lens distribution are consid¬ 
ered, the distortions introduced by strong gravitational 
lensing can be easily spotted and corrected for. For ex¬ 
ample, in the framework of the Hubble Frontier Fields 
program faint background galaxies at z ~ 7 magni¬ 
fied by up to a factor of ~ 30 by foreground massive 
galaxy clusters were dete cted down to the a bsolute mag¬ 
nitude of Muv —15.5 (jAtek et al. 11201,^ . which is al¬ 
most two orders of magnitude dim mer than the faintest 
galaxies observed in random fields (iBonwens et al.l[2ni4 
iSchenker et al. II2013I: iMcLure et al. II2013I1 . Along sim- 
ilar lines. lAlavi et al. I (I2014H used strong lensing by a 
foreground galaxy cluster to detect z 2 galaxies which 
are two orders of magnitude fainter than what is nor¬ 
mally observed at this redshift. Because strong lens¬ 
ing is so useful for the detection of faint high redshift 
sources, its statistics and properties were widely studied 
in literature (e.g., bv lWvithe et al.l (1201111 and references 
therein). However, strong lensing events are relatively 
rare with a raw probability for multi ple images at the 
highest source redshifts of o nly 0.5% (jBarkana fc Loeb 1 
120001 : iGomerford et al.l[2002ll . On the other hand, galax¬ 
ies which are not strongly lensed still undergo magnifi¬ 
cation (or de-magnification) by foreground structure and 
can be mildly magnified without being multiply imaged. 
In this case the effect of lensing can be easily overlooked 
leading to an over (or under) prediction of the number 
counts of background galaxies, and thus, to an erroneous 
estimation of their properties. The effect of the interme¬ 
diate and small magnifications on the observed properties 
of high redshift sources, such as their luminosity function, 
has not been properly addressed in the literature and is 


considered in this work for the first time. However, the 
importance of lensing events with interm ediate magni¬ 
ficati ons is being realized; for instance, (|Mason et al. I 
l2015f) mention the relevance of newly detected at z ^ 8 
sources with 1.4 < /r < 2 for the determination of the 
intrinsic luminosity function. 

Here we explore the effect of gravitational lensing by 
two types of foreground populations: (i) virialized halos 
hosting bright galaxies, and (ii) proto-clusters, i.e., non- 
virialized mildly non-linear overdensities, on the luminos¬ 
ity functions of the sub-mm and Lyman-break galaxies. 
In §2 we outline the basics of our lens model and com¬ 
pute the probability for lensing at each magnification, 
P(^), which is necessary for deriving the observed lumi¬ 
nosity functions. In §2.1 we carefully examine the case 
when the background sources are LBGs and lensing is 
done by virialized halos which host bright galaxies pay¬ 
ing particular attention to the surface brightness of the 
source and the lens galaxies. In this setup the source 
images could be hidden behind an extended lens that 
appears to be bright in the observable range of wave¬ 
lengths, which, as we find, has a crucial effect on the 
lensing statistics. In §3 we present our results, show¬ 
ing that if not accounted for, gravitational lensing with 
intermediate and small magnifications (/x < 2) can be 
responsible for errors in the derived parameters of the 
observed source luminosity function of the two popula¬ 
tions. Finally, we conclude in §4. Throughout this work, 
we adopt the sta ndard set of cosmological parameters 
(lAde et al. Il2014^ . 


2. LENSING MODEL 

A central ingredient of our calculation is the probabil¬ 
ity for lensing with each given magnification, T’(/x). We 
analytically comp ute this quantity large l y follow ing the 
approach taken bv iLima. Jain fc Devlin I (|2010bD . whose 
main assumption is that for each trajectory of a photon 
emitted by a source there is a single massive object which 
plays a dominant role in the photon’s deflection and mag¬ 
nifies the background light by a factor of /i through grav¬ 
itational lensing. In the cases when multiple images are 
expected to form we keep only the stronger magnified 
image motivated by the fact that it is easier to observe. 

Our computational approach can be summarized as 
follows. First, for a lens hosted by a halo of mass Mi 
and a source hosted by a halo of mass Ms located at red- 
shifts zi and Zs respectively the magnification at each im¬ 
pact parameter in the lens space, r, is computed. Next, 
running over all possible parameters that describe the 
lens-source system, one finds the fraction of the param¬ 
eter space, /^, which yields magnifications higher then 
/i considering only the stronger magnified image in the 
case when multiple images are expected to form. Finally, 
the lensing statistics is constructed by computing the 
probability for lensing with magnification larger than /i, 
P{> p.) = 1 — and then the probability for lensing 

with each given magnification, P(/j.) = —dP{> /x)/d/i, is 
calculated, while the total flux is conserved (< /x >= 1) 
and P(/x) is normaliz ed to unity over the entire ra nge 
of magnifications (see lLima. Jain fc Devlini (|2010bf) for 
complementary details). 
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The factor can be simply written as, 


/m = 


dzi 


d log Ms 


H{zi) 
dn 


d log Ml 


dn 


d log Ml 


N d log Ms 




( 1 ) 


where z;, Mg, Zs) is the angular cross-section for 

lensing with magnification larger than /i. If the red- 
shift distribution of sources is given, we should aver¬ 
age ffj, over Zs as well, as we do in the case of sub- 
mm gala xies taking thei r redsh ift distribution from the 
work bv iSimpson et al.~l (|2014ll . In equation ([T]) Da is 
the angular diameter distance, and H{z) is the Hub¬ 
ble constant at redshift z. The comoving number den¬ 
sity of objects, dn/dlogM and the normalization factor, 
N = J dlogMsdn/dlo gMs, are found using th e Sheth- 
Tormen mass function (jSheth fc: Tormen~III999ll . To cor¬ 
rectly estimate the number density, dn, of the objects 
with masses between M and M + dM per a logarith¬ 
mic mass interval, dlogM, one needs to specify the vari¬ 
ance of fluctuations in the matter density at each mass 
scale (which we calculate using the outputs of CAMEQ), 
and the critical overdensity at which the objects form. 
Sc- The latter quantity depends on the type of the 
objects for which we want to know the number densi¬ 
ties. As noted above, in this work we are interested in 
the effects produced by two kinds of foreground popula¬ 
tions: (1) virialized massive halos for which the critical 
overdensity is jus t the s tandard value for collapse (e.g., 
iBarkana fc Loeb I (|200IIB . and (2) non-virialized objects. 
To model the population of non-virialized objects we rely 
on the fact that the radius of a growing overdensity is 
close to its value at turnaround (the moment at which, 
in the framework of the spherical collapse model, the ex¬ 
ternal mass shell has zero velocity) during most of the 
evolution time of the overdense region. We then find the 
critical overdensity, 5c, for the objects at turnaround by 
consistently solving the spherical collapse model with rel¬ 
evant cosmological parameters at each redshift, and plug 
it into the Sheth-Tormen prescription to get the number 
density of non-virialized halos. 

Next we calculate the value of magnification at 
each impact parameter, which depends on the three- 
dimensional density profile of each lens as well as on its 
distance from the source and the observer. In treating 
objects at turnaround we consider their mean density 
in the context of the spherical top-hat collapse model. 
We follow the standard pro cedure (e.g., outlined by 
iLima, Jain fc Devlin I (l2010bll ') to find the magnification 
pattern for each overdensity. In general, the objects at 
turnaround are only mildly nonlinear having the over¬ 
densities of the order ~ 5 (with the exact value being 
redshift dependent in ACDM) which results in magnifi¬ 
cation of /X < 2. Such objects do not produce multiple 
images and thus can have only a mild effect on the ob¬ 
servable population. 

In the case of virialized halos the model is 

more complicated. _ Conve ntionally, either the 

iNavarro. Frenk fc White I (II997I1 (NFW) or the sin¬ 
gular isothermal sphere (SIS) profiles are used to model 
the three-dimensional density distribution. However, 


^ http: //camb.info 


neither of the two profiles describes lensing by realistic 
galaxies close enough. In particular, observations of flat 
rotational curves in massive galaxies show that the mass 
distri bution follows th e SIS profile within 10-20 inner 
kpc (jKochanek 1 1199411 . while the outskirts are better 
descr ibed by the NFW profile (IMandelbaum et aTTI 
1200511 . Therefore we adopt a combination of SIS and 
NFW through a piecewise form for the magnification as 
a function of the impact parameter in the lens plane 


f MS/S, r <rc 
\ fJ'NFW, r > Tc 


( 2 ) 


and connect the two profiles at a projected radius Tc 
that characterizes the extent of the luminous core of the 
galaxy (twice the half light radius). 

The last ingredient needed for equation m is the an¬ 
gular cross-section for lensing, /\Q,^{Mi, zi. Mg, Zg). To 
calculate this quantity for each lens-source pair we go 
over all r < rg, where rg is the maximal value of the im¬ 
pact parameter that allows a single lens to dominate. In 
the case when multiple images are produced we consider 
the r which corresponds to the brightest image only. By 
averaging over all the masses of lenses at each redshift, 
we first calculate the mean separation between halos, s, 
in terms of either the mean virial radius (averaged over 
halos masses at each redshift) or the typical radius of a 
structure that is turning around. In the former case we 
find s ~ 10, while in the latter case s ~ 2.5, with the 
exact value in each case being redshift dependent. Next, 
we assume that each given object dominates the lens¬ 
ing effect for the impact parameters smaller than s times 
its radius, for which we adopt either a virial radius, and 
thus rg ~ 10 Tyir, or a radius at turnaround, rg ~ 2.5 rta, 
depending on the case studied. 

The resulting P(m) is shown on the left panels of Figure 
[T]for the two source populations, i.e., LBGs at z^ = 6 and 
for SMGs at Zg = 2.6 (which represents the center of the 
current redshift sample of SMGs), and for the two lens 
populations, i.e., for virialized and non-virialized halos. 
For LBGs, a particular care should be taken when com¬ 
puting the lensing statistics since the foreground galax¬ 
ies are normally bright in the band which refers to the 
rest-frame UV of the LBGs. This effect introduces a 
suppression in the probability for lensing at strong and 
intermediate magnifications, which we also show on the 
Figure, and discus in full detail in §2.1. 

For both sets of sources, LBGs and SMGs, and in 
the case when lensing is due to the highly nonlinear 
overdense regions (i.e., virialized halos), the large mag- 
nification tail of P in) scales as oc as expected 
((Turner et al. I Il984fl , and all values of magnifications 
can be obtained, including very large // for which Ein¬ 
stein rings are produced. On the contrary, in the case 
when lensing is due to the mildly non-linear objects, P{p,) 
drops very fast with magnification, not allowing for the 
possibility of strong lensing at all. However, as we see in 
§3, the effect of the objects at turnaround on the bright 
end of the luminosity function can be comparable to lens¬ 
ing by halos when observing a sample of field galaxies 
which do not experience strong lensing. 

The effect of gravitational lensing can also be inter¬ 
preted as modification of the luminosity distances. On 
the right panels of Figure [T] we show the probability 
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Fig. 1. — Probability distribution for lensing with magnification fi (left panel) and probability distribution for 5dl (right panel) shown 
for a population of the Lyman-break galaxies (LBGs) a,t = 6 (top) and of sub-mm galaxies (SMGs) at z^ = 2.6 (bottom). For SMGs, 
the lenses under consideration are virialized halos (solid) and proto-clusters (dotted). In the case of LBGs the sources are lensed by a 
population of virialized galaxies including the reduction in P(/r) (dashed) and ignoring it (solid), as well as by proto-clusters (dotted). 


distribution for the fractional change in the luminosity 
distance due to the magnification of the flux, Sdl = 
— 1. Although the probability distribution oi 5dl 
peaks around zero (i.e., magnification has no effect in av¬ 
erage around the sky) there is dispersion in the values of 
the luminosity distances simply because along a random 
line of sight the flux can be magnified by a random fac¬ 
tor /r. This dispersion in the luminosity distance of high 
redshift objects should be carefully accounted for. In 
particular, the precision with which cosmological param¬ 
eters from a high-redshift sample of standard candles can 
be determined is expected to be affected by gravitational 
lensing. However, a more quantitative determination of 
the precision in the cosmological parameters measured 
using the high-redshift standard candles, as well as the 
effect of lensing, is beyond the scope of current paper. 

2.1. Reduced lensing probability for the LBGs behind a 
bright lens 

Even when a source is located behind a bright lens, 
it can still be separate d based on its different colors 
(jBarkana fc Loeb1l2000D . Here we apply another argu¬ 
ment, showing that if magnified enough, the source can 
outshine the foreground galaxy and be observable even 
when located behind its extended luminous core, i.e., sur¬ 
face brightness may be larger than that of the lens for 


some values of the impact parameter. This argument 
introduces a new criterion, previously ignored in the dis¬ 
cussions of the galaxies’ luminosity function, which is 
relevant in the case of the high-redshift LBGs. Since 
foreground galaxies are normally bright in the observed 
bands which correspond to the rest frame UV bands of 
the sources, some images are too faint to be seen through 
the bright part of the lens even when they are magnified. 
As a result, the probability to observe LBGs with inter¬ 
mediate and large magnifications is reduced, as we show 
below. 

In our analysis we relate to each massive halo an ex¬ 
ponential s urface brightness profi l e with a half ligh t ra¬ 
dius rc/2 (|Szomoru et al. I [20131 iKravtso-iH 120131) and 
a total UV luminosity Lj/y. The UV luminosity of 
each halo is directly proportional to the star forma¬ 
tion rate in the halo, (Luy/eig s~^ Hz~^) 7 x 

10^^ (SFR/Mq yr“^) (iMadaii. Pozzetti fc Dickinson I 
IT^. which in turn scales with its mass and redshift 
(iBehroozi. Wechsler fc Conrov I I201.H iBehroozi fc Silk I 

The lensing kernel for LBGs at Zg ^ 6—10 is dominated 
by contribution from lenses at z; ~ 1, thus implying that 
the emission wavelength of Zs 10 sources is almost an 
order of magnitude shorter than that of the lens galaxy. 
When the high-redshift sources are observed at the rest- 
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frame wavelength of 1500 A (e.g, by WFC3 camera on 
board HST), the foreground galaxy would frequently be 
observed at the wavelengths bluer than the Balmer jump. 
Ignoring lines, the spectral energy distribution (SED) 
of galaxies at the corresponding wavelength band is ex¬ 
pected to be rather flat. We use a toy model to in¬ 
clude the SED of foreground galaxies and account for 
K-corrections. Our toy model for the galaxy SED con- 
sis ts of a power-law with the observed slope according 
to iKurczvnski et al.l (j2014ll for wavelengthes within UV- 
continuu m. We mode l the B almer jump adopting obser¬ 
vations bv lOteo et al.l (j2014f l. and assume a flat spectrum 
at longer wavelengthes. A more detailed treatment of K- 
corrections in a more solid way, e.g. u sing SED models 
provided bv iBruzual fc Chariot I (|2CI03ll . goes beyond the 
scope of this paper. 

Using all the model ingredients described above we cal¬ 
culate the observed flux per unit frequency for both the 
source and the lens including the magnification of the 
source. This allows us to find impact parameters (and 
magnifications) for which the source surface brightness 
is higher than that of the lens. Such regions of the pa¬ 
rameter space contribute to the lensing statistics. In 
Figure [2] we show the observed fluxes for several lens- 
source pairs accounting also for the magnification of the 
background galaxy by the foreground galaxy. Specif¬ 
ically, we show the cases of two sources Mg = 10^^ 
Mq at Zs = 6 (left panel) and M,, = 5 x 10® Mq at 
Zs = 10 (right panel) lensed by three different lenses 
(Mi = lO^o Mq, Mi = 10“ Mq and M; = 10“ M©) 
at zi = 0.5. The masses of the deflectors were chosen 
to cover the most typical (Mi = 10^® M©) and massive 
(Mi = 10^^ ]^ 0 ) structures at that epoch. Our choice 
of Mg is motivated by the derived masses of the ha- 
los hosting Lyman-break ga laxies recently reported by 
iBarone-Nugent et al. I (j2014l l. In some cases, e.g., for an 
Mi = 10^^ M© lens at zi = 0.5 and a source of M^ = 10^^ 
M© at Zs = 6 (blue curve on the left panel of Figure [2]), 
there is a range of impact parameters for which the mag¬ 
nified source is brighter than the lens and thus can be 
observed. In this particular case, the range of impact 
parameters at which the image is visible includes r > 13 
kpc and 1.3 < r < 1.6. Note that the range of r smaller 
than the radius of the Einstein ring (here te = 1-46 kpc) 
shown with light grey curves on the Figure, is attributed 
to the weakly magnified images and is not included in the 
calculation of P{fi) plotted in Figure[T] This means that 
for r > rE the range of magnifications between /r = 8.8 
at r = 1.6 kpc and /r = 1.2 at r = 13 kpc is not observ¬ 
able and such magnifications do not contribute to our 
lensing statistics. The right panel of Figure [2] shows that 
this phenomenon is essentially the same for a different 
choice of the mass and redshift of the source. However, 
the precise range of the impact parameters within which 
the magnified source can be observed in each case de¬ 
pends on the masses and redshifts of both the lens and 
the source. 

The fact that in some cases (for some values of the 
impact parameter) the foreground galaxies outshine the 
background LBG sources affects statistics of gravita¬ 
tional lensing, as shown in Figure [T] Within the out¬ 
lined setup, probability of obtaining any magnifications 
(1 < /r < oo) is affected by this selection criterion, as can 


be seen from the Figure where F’(/r) is shown with and 
without such a reduction; however, the slope of the prob¬ 
ability distribution at large magnifications /x > 2 (strong 
lensing regime) is not strongly affected by this effect and 
scales as oc As shown in the next section, the re¬ 

duction in P(/x) affects the observed luminosity function 
of LBGs at its bright end. 


3. LUMINOSITY FUNCTION OF HIGH-REDSHIFT FIELD 
GALAXIES 

To account for the effect of gravitational lensing on the 
luminosity function, one must convolve the differential 
probability P{p) which we found in §2 with the intrinsic 
luminosity function of the population under considera¬ 
tion, 4'i„t(L). The observed luminosity function, 


'k ohs {Pobs') 





int{P — Lobs I /^)d/X, 


( 3 ) 


is modified by lensing and depending on P(/i) may have 
a completely different shape than the intrinsic one. 

The effect of lensing on the observed properties of 
sources depends on the shape of the intrinsic luminosity 
function of the population under consideration. In par¬ 
ticular, if 'ki„t(L) is flatter than L~^, then lensing has 
only a minor effect on the ob served luminosity function 
(iBlandford fc Naravan IfTO^ . Here we apply our find¬ 
ings to the LBGs at redhifts Zs > 6 and to the SMGs at 
Zg = 2.6. Although rather steep, the faint-end luminosity 
function of each one of these two populations appears to 
be flatter than the critical dependence, i.e., 'k(T) oc 
Therefore lensing is not expected to have any impact on 
the number counts of intrinsically faint sources. How¬ 
ever, at luminosities higher than a characteristic value, 
L* (or at the fluxes larger than a characteristic value 
S*), the luminosity function is expected to drop rapidly 
following either an exponential or a steep power-law. In 
this regime the effect of lensing is dramatic: it boosts 
the observed number counts of the brighter sources by 
magnifying the luminosities of intrinsically fainter ones. 
Since the abundance of intrinsically faint galaxies is much 
higher than that of the bright ones, the increase of the 
observed number counts at the bright end is striking. 
Specifically, the exponential drop of the Schechter func¬ 
tion of LBGs trans forms into a power- law with the index 
of -3 when lensed (|Wvithe et al.ll20lTI) . 

At present, the bright end of the luminosity function 
is not very well constrained for both populations that 
we co nsider here due to the substantial cosmic vari¬ 
ance (I Newman &: Davi^ I2002t ISomerville et al. 1 12004 
iTrenti &: Stiavellil 1200811 . In the particular case of the 
LBGs at Zs 9 — 10, the data from GOODS-N or 
GOODS-S WFG3/IR surveys (with the total area of ^ 
150 arcmin®) suffer from the cosmic variance of 15 — 20% 
on the overall number counts out of which the lumi- 
nous sources beyo nd L* make only a small fraction 
(|Oesch et al.l I20l4) . Similarly, in the case of Zs ~ 6 
LBGs, counted in 5 independent 20 x 7.5 GANDELS 
survey fields, a total uncertainty of 10% on the vol- 
ume density of galaxie s from cosmic variance is estimated 
(jBouwens et al.l[2014ll . In the case of the SMGs, the lu¬ 
minosity function is poorly understood at all the lumi¬ 
nosities, and its redshift evolution and S* remain uncon- 
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Fig. 2.— The top panels show the fluxes (in erg Hz“^ kpc~^ units) of foreground galaxies (dotted lines) and the fluxes of a background 
galaxy magnified by each lens (solid lines) as a function of the physical impact parameter in the lens plane. Left: Source of mass Ms = 10^^ 
Mq at redshift Zs = 6. Right: Source of mass Ms = 5 X 10 ® Mq at redshift Zs = 10 . In each case three types of foreground galaxies at 
zi = 0.5 are shown: Mi = 10 ^® Mq (red), M; = 10 ^^ Mq (green), and Mi = 10 ^® Mq (blue). The bottom panels show the corresponding 
magnification in each case in the observable part (solid lines) and non-observable part (dotted lines). Light grey lines correspond to the 
weakly magnified image in each case and are shown for completeness. 


straineifl. The first detailed surveys in the sub-mm range 
are currently on their way, for example with ALMA, and 
are expected to provide better constraints on the number 
counts of SMGs. The effect of lensing on the luminosity 
functions of both LBGs and SMGs may appear to be sig¬ 
nificant once these populations are better explored. In 
the following we analyze the effect of lensing on the lu¬ 
minosity functions for all possible magnifications, while 
paying particular attention to the impact of magnifica¬ 
tions in the intermediate range /i < 2, which is expected 
to affect most of the galaxies in the field, and comparing 
it with the case of strong lensing. 

3.1. LBGs 

The luminosity function of high-re dshift galaxies is 
commonly fitted by the Schechter form (|Schechter II1976I : 
lOesch et al.ll20l4lBouwens et al.ll2014l) 



It is characterized by a power law dependence with index 
a at the faint-end and an exponential drop at luminosi¬ 
ties higher than the critical value L*, where 'L* deter¬ 
mines the overall normalization of the number counts. 
As discussed above, the abundance of rare sources at the 
luminous end of the luminosity function, L ^ L*, is very 
sensitive to lensing due to the steepness of 'I'(L), which 
cuts-off exponentially. 

The effect of strong lensing on the luminosity function 
of the Schechter form for high-redshift so urces has been 
extens ively studied in literature (see e.g., IWvithe et al.l 
(120111) 1. It was shown that in the presence of significant 

® Although several studies provided the luminosity function of 
SMG s based on the observations by the Herschel Sp ace Observa¬ 
tory IILapi et al. Il20fll : IGonzalez-Nuevo et al. 1120121 . the quoted 
number counts at the bright-end may be overestimated as a result 
of the poor angular resolution of the satellite. Recent studies with 
ALMA, which has much higher angular resolution and is able to re¬ 
solve the SMGs at high redshifts, show that some of the previously 
unresolved bright sources are, in fact, g roups of highly clustered 
dimmer galaxies (IBussmann et al~ll2015l) . 


gravitational lensing, the luminosity function acquires a 
power-law slope with an index —3 at its bright-end, so 
that the apparent abundance of galaxies with L > L* 
is dramatically enhanced. Observing such a strongly 
distorted luminosity function would flag strong lensing. 
However, only a small part of high-redshift galaxies are 
strongly lensed, while the vast majority of galaxies expe¬ 
riences either intermediate or weak lensing. 

We start by quantifying the effect of the “common” 
lensing on the luminosity function of L BGs. Adopting 
a set of Schechter parameters found by iBouwens et al.l 
(IMl for LBGs observed at = 4 - 10 («'* = 0.44 x 
jq-o. 28 (z.- 6 ) jq -3 Mpc-3, a = -1.87 - 0.1(z^ - 6), and 
M^y = —20.97 where the UV magnitude is related 
to luminosity through M^y oc —2.5 log^^Q L*), we ap¬ 
ply our calculated P{fi) to estimate the observed lu¬ 
minosity function and bias, defined as a ratio of lensed 
to unlensed luminosity functions at a given luminosity, 

obs{L')int (L), which are shown on Figure|3]for several 
cases of maximal magnification fimax at two redshifts (6 
and 12). In particular, we show = 2, 3, and 1000, 
where the last case includes the effect of strong lensing 
with Einstein rings, including and excluding the reduc¬ 
tion in P(/i) of LBGs at intermediate magnifications dis¬ 
cuses in §2.1. Finally, we show the case of lensing by 
non-virialized halos at turn-around. 

It is clear that the reduction in P(/x) plays an impor¬ 
tant role in determining the observed luminosity func¬ 
tion of LBGs for all values of /Xmax- When measured at 
L = SOL* in the cases of Umax =2,3 and 1000, the ratios 
of the luminosity function which excludes the reduction 
in P(/x) to the one which includes it, are 49, 109, and 59 
at Zs = 6, while at Zg = 12 the numbers are 7.6, 12.5 and 
9.7, respectively. In addition, it is important to note that 
mildly lensed population, e.g., with /Xmox < 2, can still be 
described by the Schechter function with slightly differ¬ 
ent 4^* or Myy, if the latter parameter is allowed to floa10 
(where we choose to keep a fixed since the faint end num- 

® IBouwens et al.l II2014I) concluded that evolution of the criti¬ 
cal brightness Myy with redshift is insignificant and its value is 
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L/L. L/L. 

Fig. 3.— The luminosity function multiplied by L* (top panels, L*'^ is shown in units of Mpc“® and L* is calculated in erg s“^ Hz~^) 
and the magnification bias, obs{L)/'^, (bottom panels) of the LBGs at 2 :^ = 6 (left) and = 12 (right). We show the results for 
the intrinsic luminosity function (black line), and the luminosity function of LBGs lensed by a population of virialized objects including the 
reduction in P{p.) for all possible magnifications (dashed) and ignoring it (solid) for magnifications fJ, < 2 (red), p < 3 (green), fi < 1000 
(blue) as well as for LBGs lensed by non-virialized objects at turn around (black dotted). 


TABLE 1 

Summary of the errors introduced by gravitational lensing in the luminosity fungtion 
OF LBGs AT DIFFERENT REDSHIFTS. EACH ENTRY CONTAINS THE PAIR OF VALUES: (BIAS MEASURED 
AT L = 26L*, A'I'*%) INTRODUCED BY LENSING tVITH fimax = 1-3, 1.5, 1.7, 2 AND 3 DUE TO 
VIRIALIZED HALOS IN THE PARAMETERS OF THE SGHECHTER LUMINOSITY FUNCTION. We QUOTE THE 
RESULTS INCLUDING THE REDUCTION IN P{^^) (tHE UPPER SET OF NUMBERS IN EACH CASE) AND 
EXCLUDING IT (THE LOWER SET). In ADDITION, WE SHOW THE CASE OF THE LENSING BY 
PROTO-CLUSTERS WITH ALL POSSIBLE MAGNIFICATION (tHE LAST COLUMN). FITS ARE DONE FOR 
THE LUMINOSITIES IN THE RANGE OF MAGNITUDES LARGER THAN THE DETECTION LIMIT OF HST 
(for SOURCES AT Zs < 10) AND ,JWST (fOR Zs > 10) AND FOR THE BRIGHT MAGNITUDE LIMIT 
Mjjv = —24.5 CORRESPONDING TO Lmax ~ 261/*. 


Redshifts 

f-^max — 1*3 

{Lmax — 1'6 

{Lmax — 1'7 

{Lmax — 2 

{Lmax — 3 


6 

1 . 7 , 1 . 6 % 
4 . 2 , 3 . 8 % 

2 . 4 , 2.3 % 
16.6, 8.6 % 

3 . 6 , 2 . 8 % 

58.2, 13.0% 

7 . 2 , 3 . 8 % 

266.2, 18.9% 

62.1, 7.4% 

> 103, 33.7% 

1.6, 1.6% 

8 

2 . 0 , 2 . 1 % 
4 . 6 , 3 . 9 % 

3.5, 3.3% 
19.0, 9.1% 

6.4, 4.2% 
69.3, 14.0% 

15.7, 5.9% 
329.7, 20.5% 

181.2, 11.2% 

> 10^, 36.5% 

2.1, 2.7% 

10 

2 . 4 , 2 . 6 % 

4 . 8 , 4 . 0 % 

5.0, 4.5% 

21.5, 9.9% 

10.8, 5.9% 

77.9, 15.1% 

32.8, 8.4% 

401.3, 22.2% 

467.6, 15.6% 

> lO^, 39.2% 

2.8, 3.3% 

12 

2 . 8 , 3 . 0 % 
5 . 2 , 4 . 2 % 

7.1, 5.7% 
23.3, 10.4% 

17.9, 7.9% 
90.4, 16.1% 

64.5, 11.4% 
460.4, 23.3% 

> 10^, 20.8% 

> 10^, 41.9% 

4.0, 4.5% 

15 

3 . 7 , 4 . 0 % 
5 . 4 , 4 . 5 % 

12.1, 7.9% 
27.0, 11.4% 

37.1, 11.4% 
105.2, 17.5% 

166.0, 16.4% 
569.6, 25.6% 

> 10^, 30.0% 

> 103, 46.2% 

7.0, 6.0% 


her counts should not be affected by lensing). However, 
in fields which include stronger lensing events, e.g., for 
Pmax > 3, the shape of the luminosity function is signif¬ 
icantly distorted. Even for jimax = 3 the Schechter form 
develops a “secondary knee”, which (when stronger mag¬ 
nifications are included) transform s into the power law 
with the asymptotic index of —3 (jWvithe et al.l [2?)TTI) . 
Visually the effect is similar for all the redshifts at which 
the luminosity functio n has been observed, Zg ^4—10 
(|Bouwens et al.ll2014D . Lastly, the effect of lensing by 
non-virialized objects is small, but non-negligible. In 
fact, in the parts of the sky where there are no strong 
lenses it can be comparable to the lensing by virialized 
halos with pLmax 1-4 at all redshifts (as we show below 

nearly constant (in terms of the UV absolute magnitudes M^y = 
—20.97 di 0.06). However, the evolution of the faint-end slope and 
normalization with redshift was found to be significant, e.g., the 
steepening in the effective shape of the U V luminosity function was 
found to be significant at the 5.7cr level IjBouwens et al.ll2Q14lh 


in Table [ij. 

Next, we quantify how strongly lensing with fj-max < 3 
affects the fitting parameters of the Schechter function. 
When estimating the effect on the luminosity function 
it is important to set the upper limit on the luminosity, 
Lmax, which we expect to observe. This number de¬ 
pends on both the intrinsic distribution of high-redshift 
galaxies as well as on the volume of the survey. If the 
survey volume is too small, it wi ll not sample t he bri ght- 
est rarest sources. In particular. lOesch et al. I (I20I3I1 ex¬ 
pect to find sources with apparent magnitudes in the 
range 26 — 30.5 at Zg ^ 10 (which corresponds to the 
range from Muv = —17.0 to Mj/y = —21.5) using the 
Hubble Ultra-Deep Field. In this case Lmax = 1.6L* 
and the data are not expected to be sensitive to lens¬ 
ing. Fitting the Schechter function to the lensed lumi¬ 
nosity function with fimax = 1000 yields discrepancy be¬ 
low 1% in the value of 4** when compared to the intrinsic 
case and keeping a and L* fixed. The effect of lensing 
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on the parameters of the luminosity function starts to 
manifest itself when the brightest observed sources have 
Mjjv < —22.5. In this regime, by comparing the strong 
lensing case to the no-lensing case we get a ^ 2% dis¬ 
crepancy for LBGs at Zg = 15 in the value of di'* (with 
even a weaker effect on LBGs at lower redshifts), whereas 
the effect on the critical luminosity L* is still negligi¬ 
ble. When the brightest sources reach Muv ~ —23.5, 
in which case L^ax lOL*, strong lensing starts to af¬ 
fect the shape of the luminosity function and the errors 
in fitting parameters are /\^*{nraax = 1000) ^ 6.0%, 
= 5) ~ 2.3% and /^'l>*{pLmax = 3) - 2.0% 
for the source population at Zg = 15 when the reduction 
in P{p) is taken into account (and 30%, 6.2% and 4.2% 
respectively when it is ignored); when the brightest ob¬ 
jects are of Muv ~ —24 {Lmax ~ 16L*) the errors reach 
= 1000) - 30%, I^'i>*{lJimax = 5) ~ 11% and 
/M^*{lirnax = 3) ~ 7% with the reduction in P{^JL) (and 
44%, 21.6% and 14.3% when it is ignored). In the fol¬ 
lowing we adopt the value Muv = —24.5 for the bright¬ 
est observable galaxy {Lmax/L* ^ 26) to put the upper 
limit constrains on the distortion of the luminosity func¬ 
tion of LBGs due to lensing, and summarize the results 
in Table m For this limiting value of Muv we estimate 
the discrepancy in the normalization of the Schechter 
luminosity function (the discrepancy in M^y is negligi¬ 
ble reaching ~ 2% for Zg = 15 with ^max = 3) at a 
range of redshifts relevant for the observations with HST 
and JWST for ^.max = 1-3, 1.5, 1.7, 2 and 3 including 
and excluding the reduction in P{p). For higher val¬ 
ues of Umax, he., when observing a field which includes 
strongly lensed (for example, the case of fimax = 1000 
in Figure [3]) bright sources {Lmax ^ 26L* and higher), 
the effect of lensing on the luminosity function is ap¬ 
parent and the resulting luminosity function cannot be 
fitted by the Schechter form. Our predictions for the 
JWST redshift range {zg = 12 and 15) are based on the 
extrapolation of the results for the intrinsic luminosity 
function. In addition, we quote the results for the lens¬ 
ing with proto-clusters for which the discrepancy in the 
Schechter parameters appears to be comparable to the 
case of lensing by virialized halos with 1.3 < Umax <1.5 
at = 6 — 15. In agreement with Figure [31 the results 
in Table |T] show that the effect of lensing is manifested 
stronger in the case when there is no reduction in P(/i) 
due to the surface brightness argument. For instance, 
for sources at z = 8 the error in the normalization of the 
luminosity function drops from ~ 21% in the case when 
the surface brightness arguments are ignored to ~ 6% 
when the reduction is accounted for. Therefore, the non- 
realistic case with no reduction in P(/r) can be used to 
set an upper limit on the effect of lensing. 

3.2. SMGs 

Next, we address the effects of lensing on the luminos¬ 
ity function of SMGs in the redshift range z = 1 — 4. 
Wide field surveys conducted by the Herschel Space Ob¬ 
servatory (Pilbratt et al. 2010) and the South Pole Tele¬ 
scope (Garlstrom et al. 2011) have been very effective at 
discovering gravitationally lensed galaxies in large num¬ 
bers (e.g., Negrello et al. 2010, Vieira et al. 2013). Lens 
models based on high resolution imaging with the Sub¬ 
millimeter Array (SMA) (Bussmann et al. 2012, Buss- 
mann et al. 2013) and ALMA (Hezaveh et al. 2013) are 


now becoming available for a substantial portion of these 
objects. Therefore, it is timely to apply our approach to 
this population of galaxies and quantify the effect that 
gravitational lensing may have on its luminosity function. 

As mentioned above, the luminosity function of the 
SMGs is not very well constrained at the moment. Thus 
there exist a variety of functions th at can equally well fi t 
the number counts observed so far ([Karim et al. 11201^ . 
The two most popular fits which we use here are the 
Schechter function and a broken power law. First, the 
Schechter form is generically the same one as is used in 
the case of the LBGs, 


dn 

dS 


N* 



—a 

exp {-S/S *), 


( 5 ) 


with S being the observed flux, N* the normalization, 
S* the critical flux at which the number counts start to 
decline exponentially, and a the faint end slope. Another 
frequently used functional form is the broken power law 

dn i N*(4ry^\ for S < S* 

= to, s>s- 

with TV* being the normalization. S'* the characteris¬ 
tic flux and /3i and /32 the slopes of the two power- 
la ws. In Figure IS we show the number counts derived 
bv iKarim et al. ll2013ll using the counts of faint sub-mm 
galaxies in the 870-/rm band of ALMA together with sev¬ 
eral fits to the luminosity function. In particular we show 
three cases for the Schechter function: 


1. Best fit model used bv iKarim et al. I (1201311 with 
S* = 8 mJy, N* = 424 deg“^ and a = 1.1, 

2. Flat fit with S* = 10 mJy, N* = 600 deg“^ and 

a = 0 . 1 , 

3. Steep fit with S* = 7 mJy, N* = 600 deg“^ and 
a = 1.9, 

and three cases for the broken power law fit: 

1. Fit 1, N* = 15 mJy“^ deg“^, S* = 8 mJy, and 
Pi =2 and P 2 = 5, 

2. Fit 2, iV* = 20 mJy“^ deg“^, S* = 8 mJy, and 
Pi =2 and P 2 = 6.9, 

3. Fit 3, N* = 25 mJy“^ deg“^, S* = 8 mJy, and 
Pi =2 and P 2 = 18. 

As we see form Figure |3| where we plot the intrinsic 
and the lensed with ^max = 1000 luminosity functions 
of SMGs together with the measured number counts 
(IKarim et al. ll2M^ . the effect of lensing starts to man¬ 
ifest itself at 5" = S'* in the case of the broken power law 
fits, while the Schechter fits do not show deviation until 
S ~ IQS'*. The strength of the effect introduced by lens¬ 
ing depends on the steepness of the luminosity function 
at its bright-end. For instance, in the case of the power- 
law dependence, for our steepest broken power-law fit 
(Fit 3) the bias reaches the value of 10 at S ~ 2S*, 
while for the flattest Fit 1 it reaches same value only at 
S' ~ bos'* . 

In Figure [4] (bottom row) we compare the intermedi¬ 
ate and strong lensing of SMGs for the steepest models 
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SMGs, Schechter LF 



SMGs, Power-Law LF 



SMGs, Power-Law LF 
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Fig. 4. — Top row: The intrinsic (dotted) and lensed (solid, fimax = 1000) luminosity function of the SMGs number counts (in mJy ^ 
deg“^ units) shown for the Schech ter (left) and the bro ken power law (right) fits with the color code marked on each panel. The black data 
points are taken from the paper bv IKarim et al. ]120T3). Note that the number counts at the highest measured flux value, S ~ 10 mJy, can 
be underestimated due to the effect of cosmic variance. Therefore, we do not strictly account for this point when fitting the data. Bottom 
row: The luminosity function (top panels) and the magnification bias (bottom panels) of the SMGs fitted by our steepest fits. We show 
the intrinsic luminosity function (black line), and the luminosity functions lensed by a population of virialized objects for magnifications 
11 <2 (red), fi < 3 (green), fi < 1000 (blue) as well as by proto-clusters (black dotted). 


TABLE 2 

Summary of the errors introduced by gravitational lensing in the luminosity function of SMGs at 
Zs = 2.6. Each entry contains the values of (Bias measured at 5 = 25S*, AAi*%, AS*%) introduced by 
lensing due to virialized halos with flmax = 1.3, 1.7, 2 AND 3, AS WELL AS PROTO-CLUSTERS, IN THE 
PARAMETERS OF THE SCHECHTER LUMINOSITY FUNCTION FOR THE INTRINSIC FUNCTIONS SHOWN IN EiGURE ThE 
FITS ARE DONE FOR THE LUMINOSITIES IN THE RANGE 3 mJy < S < 255*. 


Model 

f-i'max — 1*3 

f-^max — 

l-^max — 2 

fJ'tnax — 3 

. ta 
^all 

Flat 

4.9, 2.7%, 5.9% 

20.9, 8.2%, 9.9% 

56.4, 11.9%, 12.7% 

511.9, 20.2%, 20.0 % 

1.1, 0.5%, 0.2% 

Best fit 

5.5, 5.2%, 6.3% 

30.4, 15.8%, 11.3% 

96.6, 23.1%, 15.0% 

> 10^, 38.5%, 25.6% 

1.1, 0.7%, 0.2% 

Steep 

6.1, 7.5%, 6.6% 

41.5, 23.0%, 12.5% 

151.3, 33.2%, 17.1% 

> 10^, 53.3%, 30.7%, 

1.1, 0.9%, 0.3% 


in which the effect of lensing is the strongest (i.e., our 
steep fit in the case of the Schechter function and Fit 3 
for the broken power law), showing the impact of lens¬ 
ing with ^max = 2, 3, and 1000 on the number counts 
of SMGs in the case when lenses are the virialized ha¬ 
los as well as the case of lensing by proto-clusters. We 
also show the bias, defined as the ratio of the observed 
to the intrinsic number counts, {dn/dS)^^^ / {dn/dS)^^^. 
In addition, the effect on the derived parameters of the 
SMG Schechter luminosity function and bias are further 


elaborated in Table [2] for the cases of lensing by halos 
with ^rnax = 1-3, 1.5, 1.7, 2 and 3 as well as lensing 
by proto-clusters. As in the case of LBGs, the effect 
on the parameters of the observed luminosity function 
depends on the upper limit of the luminosity that we 
expect to measure. Here we adopt the upper value of 
Smax = 255*. Both the Figure and the Table show that 
in the case of SMGs at Zg =2.6, lensing by proto-clusters 
has a negligible effect due to the relatively low redshift 
of the sources. In particular, the bias in the case of the 
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sources with a power-law form of the intrinsic luminos¬ 
ity function lensed by proto-clusters (bottom right panel 
of Figure 2]) reaches 1.04 at the bright-end, while in the 
case of viriallized halos with ^.max = 2 it is ^ 33. By 
comparing the numbers in Table [21 we find that the ef¬ 
fect of proto-clusters in this case is negligible even when 
compared to lensing by halos with jSrnax = 1-3, unlike in 
the case of the LBGs. 

When the SMGs intrinsic luminosity function is given 
by the broken power law, its shape is completely dis¬ 
torted even after it goes through a relatively mild lensig 
(m < 3), as we see in Figure [4] where the effect on Fit 
3 is shown. (However, the intermediate lensing has al¬ 
most no effect when a flatter power law, e.g., our Fit 1 
or Fit 2 with ^2 = 5 and 6.9, respectively, is considered). 
In this case a different type of dependence (e.g., a three 
power-law piece-wise function) would be needed to fit 
the data. More specifically, for lensing with some ^imax, 
a second “knee” develops around Sobs ~ Mmax x S*. At 
the flux values S* < S < Umax x S* the luminosity func¬ 
tion tends to develop slope equivalent to the one in the 
case of strong lensing (which is —2 in the case when num¬ 
ber counts per unit flux are concerned), while at fluxes 
higher than fXmax x S* the luminosity function retains 
its intrinsic power-law dependence. Thus, the effect of 
lensing can be easily corrected for in this case. 

4. CONCLUSIONS 

We considered the effects of gravitational lensing on 
the luminosity functions of Lyman-break galaxies at 
z ~ 6 — 10 and sub-mm galaxies at z ^ 2 — 4 in 
the regime when a single lens has a dominant contribu¬ 
tion. In particular, we made predictions for the galaxies 
“in the field”, i.e., in surveyed volumes which are not 
strongly lensed, focusing on the intermediate range of 
magnifications (/x < 2) produced by two types of lenses, 
namely virialized halos and non-virialized proto-clusters, 
and comparing the results to the cases which do include 
strong lensing. The regime of intermediate lensing was 
not studied in the literature and may appear to be impor¬ 
tant for future precision analysis of these two populations 
of sources when observed by instruments such as HST, 
JWST, ALMA and Herschel. 

We find that when a population of Lyman-break galax¬ 
ies with a Schechter intrinsic luminosity function is sub¬ 
ject to lensing with magnifications < 3, the ob¬ 
served luminosity function can still be fitted by the 
Schechter form, although with modified Schechter pa¬ 
rameters. Therefore, errors in derived intrinsic parame¬ 
ters are expected if de-lensing of the field galaxies is not 
done carefully enough. It appears however, that the com¬ 
mon effect of lensing with intermediate magnifications is 
negligible when existing surveys with the field of view 
of HST are considered because of the significant cosmic 
variance (~ 10% at z 6 and ^ 20% at z ^ 10). In par¬ 
ticular, we found that the overall effect of lensing (includ¬ 
ing the stro ng lensing effects) on the current Schechter 
parameters (iBouwens et al.ll2014f) is well below 1% when 
only sources with UV magnitudes above Muy —21.5 
are present (as is expected for the most luminous objec ts 
in the HST surveys at z ^ 9 — 12 (jOesch et al. 11201311 1. 
However, the effect of lensing starts to be increasingly im¬ 
portant if intrinsically brighter sources are found in the 
sample, which may be the case in the future with wider 


surveyed areas. In this regime, the errors introduced 
by the weak and intermediate gravitational lensing may 
become significant, e.g., lensing with maximal magnifica¬ 
tion of 2 of a population of bright sources at z = 12 may 
result in more than 20% discrepancy (and more than 40% 
discrepancy when the magnifications in the range /x ^ 3 
are considered) in the normalization of the luminosity 
function if sources with Muy < —24.5 are present (while 
the effect on critical luminosity remains at a percent level 
and we keep the faint-end slope fixed). If indeed intrin¬ 
sically bright sources with Muy < —24.5 exist so early 
on in the cosmic history, they would be valuable probes 
of the early Universe. These rare objects are expected to 
form in highly overdense regions, to produce stars very 
early, and to start heating and ionizing gas earlier than 
in average over the Universe. Consequently, during the 
epoch of reionization (z > 7) such regions would cre¬ 
ate big bubbles of ionized gas over the otherwise neutral 
background, thus imprinting strong signature in the red- 
shifted 21-cm signal of neutral hydrogen which can be 
probed through tomography by radio telescopes such as 
the Square Kilom eter Array (|Loeb fc Furlanett^ 120131 : 
lAhn et al. Il20l4 l. 

An important point that we revealed in this paper is 
that particular care should be taken in the case in which 
Lyman-break galaxies are lensed by massive halos that 
themselves host bright galaxies. In this case it is crucial 
to account for the magnification of the surface brightness 
of the source and compare it to that of the lens when de¬ 
ciding which sources are, in fact, observable and comput¬ 
ing the lensing statistics. We showed that the range of 
impact parameters (and thus magnifications) for which 
the high-redshift sources are observable is reduced when 
the constraints from surface brightness profiles are ac¬ 
counted for. This leads to a suppressed probability for 
obtaining intermediate and strong magnifications which 
could result in up to a two orders of magnitude deficit in 
the luminosity function at high redshifts and high lumi¬ 
nosities. 

Comparing lensing by virialized halos, which are 
strongly non-linear systems, to the role of proto-clusters, 
which we assumed to be overdense regions at the point 
of turnaround, we found that for a population of Lyman- 
break galaxies at z = 6 — 15 the effect of proto-clusters is 
similar to that of lensing by massive halos with the cutoff 
1.3 < /imaa; < 1-5 in magnification. For example, proto¬ 
clusters can introduce ^ 5% discrepancy in the normal¬ 
ization of the number counts of sources at z = 12, while 
lensing by virialized halos with fj^max = 1-3 contributes 
only ^ 3% when the suppression in lensing probability 
due to the surface brightness is taken into account (and 
~ 4% when it is ignored), and lensing with fimax = 1-5 
contributes ^ 6% (~ 10%). 

The effect of lensing on the luminosity function of sub- 
mm galaxies at z 2.6 was also considered in this work. 
Since the intrinsic luminosity function of this popula¬ 
tion is highly unconstrained at the moment, the effect of 
lensing could appear to be either significant or negligible 
once this population is better explored, depending on the 
steepness of the luminosity function at the bright-end. In 
this work we considered two forms of the intrinsic lumi¬ 
nosity function for this population, i.e., the Schechter 
fit and the broken power law dependence. In the lat¬ 
ter case, lensing with intermediate magnifications may 
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cause significant distortions of the shape of the luminos¬ 
ity function, and thus the effect of lensing can be easily 
corrected for. The former case is qualitatively very simi¬ 
lar to that of Lyman-break galaxies, with the sole differ¬ 
ence that there is no reduction in the lensing probability 
due to the large surface brightness of an extended lens. 
This is because the foreground lenses, which are drawn 
from common galaxies, are not expected to be bright 
at sub-millimeter wavelengths. The discrepancy in the 
Schechter parameters introduced by lensing by virialized 
halos with intermediate magnifications (/i < 2) is ^ 23% 
in the normalization and ^ 15% in the critical luminos¬ 
ity when the best fit to current data is considered at flux 
values 3 — 200 mJy. Finally, we found that the effect of 
lensing by proto-clusters on the galaxies at z ~ 2 — 4 
is negligible when compared to the lensing with inter¬ 
mediate magnifications by virialized halos. For exam¬ 


ple, for the best-fit intrinsic parameters that we consid¬ 
ered above, proto-clusters introduce ~ 1% discrepancy in 
the normalization of the number counts, while lensing by 
virialized halos with fimax = 1-3 contributes more than 
5%. 

Lastly, we point out that magnification can have an 
important effect on the determination of the luminosity 
distances to the high-redshift standard candles. As a 
result of randomly distributed foreground structure the 
luminosity distance to each source is subject to an er¬ 
ror, which, if unaccounted for, will bias the values of the 
deduced cosmological parameters. 
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